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I. INTRODUCTION 

Recently, there has been a substantial effort to study the properties of cold and warm 
quark matter. Such studies are important to understand the properties of the astrophysical 
compact objects like neutron stars and pulsars. There has been lot of experimental efforts 
like Einstein laboratory, ROSAT, CHANDRA and XMM, where various measurements are 
performed to understand the properties of neutron stars l|, |2|. 

There is a possibility that at the core of neutron stars the density may go upto 5 ~ 6 times 
the normal nuclear matter density where the matter is not expected to be in the hadronic 
phase. In fact, under such a scenario one expects that it would be more appropriate to 
describe the core of such dense stars as degenerate quark matter [l| which is our main 
interest in the present work. 

It is known that the newly born neutron stars cool via the emissions of neutrinos and 
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antineutrinos within few minutes involving neutron beta decay and its inverse reaction [3|-|5|. 
This coohng behavior is related to the transport properties of the medium or thermodynam- 
ical quantities like specific heat, entropy etc. All these quantities are in turn related to the 
neutrino scattering rate or mean free path (MFP). 

nn 

Such emissions can take place either by the direct [3H5[ or by the modified URCA processes 
{4, 6|. The direct URCA reactions can proceed in neutron rich matter if the ratio of the 
proton number density to the total baryon number density exceeds a critical value which 
follows from the energy and momentum conservation properties 3|- In modified URCA 
process another nucleon catalyzes the reaction to occur under situations where the direct 
URCA reaction is forbidden. 

For quark matter, similar emission processes proceed through the decay of d quark or the 
scattering of u-d quarks. Analogously these reactions are named as 'quark direct URCA' 
processes which have been studied in detail by Iwamoto 4]. Our main focus here is to cal- 
culate the neutrino MFP, as mentioned earlier, for cold and warm QCD matter. Previously 
MFP of quark matter was derived in Ref.js] where the calculations were restricted to the 
leading order by assuming free Fermi gas interactions. Here we plan to go beyond leading 
order to include 'plasma' or 'quasi particle' effects arising out of the interacting ground state. 
This, as we show, gives rise to logarithmic corrections to the MFP in contrast to what one 
expects from the usual Fermi liquid theory (FLT) js, ol. Similar non- Fermi liquid behavior 
of various thermodynamical quantities have recently drawn significant attention. 

For example, in (l|,llO|, the authors have computed the leading contribution to the interac- 
tion part of specific heat and entropy when temperature is much smaller than the chemical 
potential of quark matter. It has been shown, in case of specific heat, instead of (9(T^) 
term, an anomalous T In T term appear as corrections to the leading order which becomes 
dominant for temperature smaller than the Debye mass. This type of nonanalytic terms are 



known as non-Fermi liquid (NFL) corrections 



10|-ll2|. Ref.lnl examined NFL effects 



in the normal phase of high density QCD matter both using the Dyson-Schwinger equation 



and the renormalization group theory. In [12| magnetic susceptibility has also been shown 
to receive correction (9(T^lnT). Being motivated by these series of works we undertake 
the present investigation to estimate and see the consequences of such effects in the case of 
neutrino MFP. 

In FLT, quarks are treated as quasiparticles and their energy (E) is regarded as a func- 



tional of the distribution function 



BI9, 



13Nl6|. FLT is restricted to the low- lying excitations 



near the Fermi surface, where the lifetime of quasiparticles are long enough. Therefore, it 
is an important tool to study the properties of nuclear (or quark) matter. In Ref . jl7 1 . it 
has been argued, that the exchange of dynamically screened transverse gluons introduces 
infrared divergences in the quark self-energies that lead to the breakdown of the Fermi liquid 
description of cold and dense QCD in perturbation theory. A detailed study of non-Fermi 
liquid aspects of the normal state was presented in Ref. {2!. There the spectral density, 
dispersion relation and width of quasiparticles with momenta near the Fermi surface were 
derived at T = by implementing a renormalization group resummation of the leading 
logarithmic infrared divergences associated with the emission of soft dynamically screened 
transverse gluons l|, l2|. 

We have already mentioned that such anomalous corrections are ultimately connected to 
the absence of magnetic screening of gluons via. Landau damping. One such calculations 
which we find to the most relevant for the present purpose was performed in js]. It was 
shown that the emissivity receives logarithmic corrections is enhanced due to non-Fermi 
liquid effects. It might not be out of context here to recall two important works on neutrino 
mean free path in QED plasma. One is due to Tubbs and Schramm and the other is 



done by Lamb and Pethick 19|. In 18[, the resultant mean free path was calculated in the 
neutranized core and just outside the core. It is concluded in 19j, that neutrino degeneracy 
reduces the neutrino mean free path which suggest that neutrino may flow out of the core 
rather slowly. 

In our work we show that the neutrino mean free path receive logarithmic corrections 



where t 
the 



le dressed gluon propagator is used instead of bare propagator. Infact, we generalize 



19| results by incorporating the NFL corrections for the quark matter. The 



corrections to the MFP for degenerate and non-degenerate neutrinos, as we shall see, will 
involve different powers of as- 

It is well known that in the interior of the star the neutrino absorption and neutrino 
scattering can take place. These two types of mean free paths are related to different 
physical phenomenon. Absorption mean free path gives the transport properties, which 
characterizes the rate of momentum transport. This also enters into the expression for the 
diffusion coefficient. While, the scattering mean free paths are related to the relaxation 
time, which characterizes the rate of change of the neutrino distribution function {4]. To 
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obtain the total mean free path one can write [2 



1 1 _J_ , . 

7. s rait \ ) 



Itotal jabs Jscatt 

mean mean mean 



II. MEAN FREE PATH 



In this section we calculate neutrino mean free paths in quark matter including NFL 
corrections both for degenerate and nondegenerate neutrinos. When the neutrino chemical 
potential (fiy) is considered to be much larger than the temperature (T), the neutrinos 
become degenerate and for nondegenerate neutrinos fi^, <^ T. In our model the Lagrangian 
density is described by jst 

CwAx) = ^l,{x)J^{x) + H.C. (2) 

where the weak coupling constant is G ^ 1.166 x 10"^^ in MeV units, and and J'^r are 
the lepton and hadron charged weak currents, respectively. The weak currents are 

l^{x) = 67^(1 -75)z/e + /i7;,(l -75)1^^ + ...., (3) 
Jl^{x) = cose,uY{l-^5)d + sme,uYil-l5)s + (4) 



where 6c is the Cabibbo angle (cos^ 6'c — 0.948) [21 1. 

The mean free path is determined by the quark neutrino interaction in dense quark matter 
via. weak processes. We consider the simplest /3 decay reactions; the absorption process 

c? + z/g — )■ ti + e~ (5) 

and other is its inverse relation 

u + e~ d + (6) 

The neutrino mean free path is related to the total interaction rate due to neutrino 
emission averaged over the initial quark spins and summed over the final state phase space 



and spins. It is given by j4| 

1 g f d^pd 1 f d^pu 1 f d^pe 1 



f <£p^J_ f (£p^J_ f £P^J_f^.A.4fp , p _ p _ 

XE^,T) 2eJ (27r)32Ej (27r)3 2i?J (27r)3 2Ee^ ^ ^ 

\M\^{n{pd)[l - n{p^)][l - n{pe)] - n{pu)n{pe)[l - n{pd)]}, (7) 
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where, g is the spin and color degeneracy, which in the present case is considered to be 6. 
Here, E, p and Up are the energy, momentum and distribution function for the corresponding 
particle. |Mp is the squared invariant amplitude averaged over initial d quark spin and 



summed over final spins of u quark and electron as given by 



(8) 



Here, we work with the two flavor system as the interaction involving strange quark is 



Cabibbo suppressed 



a. 



A. Degenerate Neutrinos 

We now consider the case of degenerate neutrinos i.e. when /i^ ^ T or in other words we 
consider trapped neutrino matter. So in this case both the direct (Eqj5]) and inverse (Eqj6]) 



processes can occur and both the terms in Eq.Q under curly brackets {4] are retained. 
Consequently, the [3 equilibrium condition becomes + fJ'u = /^m + /^e- Neglecting the 
quark-quark interactions and by using Eqs.([7]) and ([S]), for the mean free path one obtains 

= -^G^cos^6'c / d^pd I d^pu j d^peil - cos 9du)il - cos 9ue) 

X 5\Pd + P,-Pu- Pe)[l + e-^^''^-^^^]n{pd)[l - n{p^)][l - n{p,)] (9) 

In the square bracket, the second term ^-Pi^i'-t^i') jg flue to the inverse process (Eqj6]). Since 
the masses of u, d quark and electron are very small, one can neglect the mass effect on the 
mean free path. To carry out the momentum integration we define p = \pd + Pu\ = \Pu + Pe\ 



as a variable. Following the procedure, described by Iwamoto {4] one has 

pdp 

sinOdudOdu = — rr. — (10) 
p'^ — 2p^p'l + pi 

(1-C0SM(1-C0SM ^ , ,n , , , , , , 11 

^PfKd)pAv)Pfiu)pf{e) 

and 

d^Pd = 2tx sine dude duP^dpd 

Pf{d) dpd 
= 2'K—--pdp-—dEd 
Pf{u) dEd 

= 2nP-I^pdp^du (12) 
Pfyiy) duj 

d'p^ = 2nPl^^^dE,^du: (13) 
p doj 
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where we denote the single particle energy E^^u) as u. For the free case dp/du is the 
inverse quark velocity. It is well known that this slope of the dispersion relation changes in 
matter due to scattering from the Fermi surface and excitation of the Dirac vacuum. The 
modified dispersion relation can be obtained by computing the on-shell one-loop self-energy. 
For quasiparticles with momenta close to the Fermi momentum, the one-loop self-energy is 



dominated by the soft gluon exchanges 



22 



23| 



22| . The quasiparticle energy cu satisfies the relation 



Ep{uj) + ReS(ci;,p(ci;)) 



(14) 



where we have approximated only to the real part of sel 



■ energy, since the imaginary part 



of S turns out to 
be obtained in 



De negligible compared to its real part 



23 
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24| 



12| . The detailed analysis can 



In the relativistic case, when the Fermi velocity f/ is close to the velocity of light c, 
exchange of magnetic gluons become important. In the non-relativistic case it is suppressed 
by a factor {v /cY with respect to exchange of electric gauge bosons and is usually neglected. 
The magnetic interaction, as stated before, is screened only dynamically and the problem 



can be written as 



10- 



13, 



remains for the static gluons [22 
tion has to be performed 



22| 



25| . Therefore to obtain a finite result, a suitable resumma- 



. The analytical expressions for one-loop quark self-energy 



127r2 



iu - /i) In ( — ^ ) + i 

\u — jJiJ 



12-K 



\UJ — /i| 



(15) 



It exhibits a logarithmic singularity close to the Fermi surface i.e. when w — ?■ /i. Thus the 
ong ranged character of the magnetic interactions spoils the normal Fermi-liquid behavior 



y, Il3l-ll6|. The breakdown of the Fermi liquid picture is associated with the vanishing 
of the discontinuity of the distribution function at the Fermi surface {2!, [l7|. This non- 
perturbative nature of the self energy gives rise to the non- Fermi liquid behavior. Here, niD 
is a cut-off factor and should be an order of the Debye mass. Differentiating Eq. ffT^ with 
respect to p, we obtain at leading order in — is 



dp{uj) 
dui 



1+^inr^^ 



(16) 
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where is the strong couphng constant, Cp = {N^ — l)/(2A'^c) and N^, is the color factor. 
Using Eqs. fll6p . (jH]) and flimi3p . the neutrino mean free path can be determined for two 
conditions. For \pf{u) —pf{e)\ > \pf{d) — pf{iy)\ 

1 



abs,D 
mean 



l + ^lnf^^ 



\ T J 



(17) 



To derive above Eq.f |T7|l . we use the result of integral 

POO poo POO 

/ dEd / dE^ / dE,[l + e-^(^''-^'^V(P<i)[l - ^iPuW - niPe)]6{Ed + E, - E^ - E,) 
Jo Jo Jo 



Similarly, for \pf{d) — pf{iy)\ > \pf{u) — pf{e)\, the corresponding expression for mean 
free path can be obtained by replacing ^ f^d and /ig O //^ in Eq. lITTI) . Since quark and 
electron are assumed to be massless, the chemical equilibrium condition gives Pf{u)+Pf{e) = 
Pf{d) +pf{h'), which we use to derive Eq.l fTTl) . 

The other major contribution to the mean free path arises from quark-neutrino scattering. 
The neutrino scattering process from degenerate quarks are given by 



(19) 



for each quark component of flavor i{= u or d). The scattering mean free path of the 
neutrinos in degenerate case can be calculated similarly as evaluated by Lamb and Pethick 



m 



19| for electron- neutrino scattering. Assuming rriqjpf^ <^ 1 and including the non- Fermi 



liquid correction through phase space, the mean free path is given by 



1 



scatt,D 



I 



16 



ng^ao 



(E,-/i,)2 + 7r2T2 





[l + ^^"^ In 






I Sn 


K T )\ 



A(x, 



(20) 



and is the vector and axial vector coupling constant 
^. In Eq.f l20p if we drop the color factor and the second 



Here, ruq- is the quark mass. C 
was given in TABLE-II of Ref. 
square bracketed term, we obtain results reported in 19| for dense and cold QED plasma 
with rriq replaced by nie. In Eq. fl20|) the constants cxo = ^G^mgjn [l8| and Ug- is the number 
density of quark, is given by 

d'^p 1 



n„ 



(27r)3e^(-^9«"^9.) 



(21) 



where 2 is the quark spin degeneracy factor. Exphcit form of A{xi) can be written as 



4 Min(/it.,/igJ 

3 f^qi 



{C'+Ci){2 + -xi+2Cv^CA.x 



and Xi = ^iu/Hq^ if ^-u < /^q, and Xi = figj^i^ if fiu > fJ'q, 



B. Nondegenerate Neutrinos 



HQ 



(22) 



We also derive mean free path for nondegenerate neutrinos i. e. when ^„ <^ T. For nonde- 
generate neutrinos the inverse process (jH]) is dropped. Hence we neglect the second term in 
the curly braces of Eq.(I7]). In this case, only those fermions whose momenta lie close to their 
respective Fermi surfaces can take part in a reaction. It is to be mentioned here, if quarks 



are treated as free, as discussed in H 



29 



30|, the matrix element vanishes, since u, d quarks 



and electrons are collinear in momenta. The inclusion of strong interactions between quarks 
relaxes these kinematic restrictions resulting in a nonvanishing squared matrix amplitude. 
Since the neutrinos are produced thermally, we neglect the neutrino momentum in energy- 
momentum conservation relation {4]. This is not the case for degenerate neutrinos where 
Py^T and therefore such approximation is not valid there. By doing angular average over 
the direction of the outgoing neutrino, from Eq.(|8]) the squared matrix element is given by 

B 



|M|2 = COS^ e,p){Vd ■ Pu){Vu ■ Pe) 

= 64G2^os2^cPj^i./ie^"'^' 



TT 



(23) 



where V = {l,Vf) is the four velocity. To calculate Eq. fl23|) we have used the chemical 
equilibrium condition fid = fJ'u + f^e and also the relations derived from Fermi liquid theory 
are given by^ 



vf = 1 — 

Putting |Mp in Eq.Q we have, 
1 



271 



Cpas 



IT 



(24) 



,ND 



5\Pd + Pu- Pu - Pe)n{pd)[l - n{p^,)][l - n{p,)] 



(25) 
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Neglecting the neutrino momentum in the neutrino momentum conserving 5 function, the 
integrals can be decoupled into two parts. Following the procedure of Iwamoto 4], the 
angular integral is given by 

,2 



A 



dQd / dQu / dQeS {pd - Pu - Pe 



(26) 



and the other part 



B 



>dpd 



■dEd 



2 dpu j„ 



pldE, 



6{Ed + E,-Eu- Ee)n{pd)[l - n{pu)][l - n{p,)] 



(27) 



Changing the variables to Xd = {Ed — Hd)f^, Xu = —{E^ — jJiu)P and Xe = —{E^ — /ie)/3 
and denoting the single particle energy E^i^d) as u we have from Eg. ( 1271) 

dpd{uj) dpu{uj) 2 2 2; 



B 



duj du 



-PdPuPe^i^d + Xu + Xe + [3 E^)n{x d)n{-Xu)n{-Xe) (28) 



As the contribution dominates near the Fermi surfaces, extension of lower limit is a reason- 



able approximation 28|, l31 1 . 



Using Eq. f|T6|l and perform the integration of Eq. fl28l) following the procedure defined in 



12 



28 



3l|, we have 



1 + ^lnf^) 
Sn \ T J 



2(1 + e-/5^' 

Using Eqs. (!25l) . (!26|) and (!29|l . the mean free path at leading order in T//i is given by 
1 



(29) 



abs,ND 



-G cos 9c fid l^u Aie TT^V — ^FT ^ + o 



(30) 



The first term is known from j4| and the additional terms are higher order corrections to 
the previous results derived in the present work. 

For the scattering of nondegenerate neutrinos in quark matter, the expression of mean free 
path was given by Iwamotoj^. We incorporate the anomalous effect which enters through 
phase space modification giving rise to 



iscatt,ND on ^ \ m I \ ii 

Imean ■^^ ^'"'qi^ ^ H'l 

Here, we have assumed mg-fpf- ^ 1 and the constants ctq and number density rig. defined 
earlier. 



20 



(E,\^(E, 



CpOs 



In 



T ) 



(31) 
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TABLE I: Comparison of neutrino mean free path in quark matter for the case of degenerate and 
nondegenerate neutrinos. Mean free path given in unit of meter, without NFL correction and 
including NFL correction for the case of absorption and scattering processes. 





Absorption 


Scattering 




Without Correction 


With Correction 


Without Correction 


With Correction 


Degenerate 


3.62 


1.85 


15.88 


8.13 


Nondegenerate 


7.68 


3.93 


1585.18 


811.06 



III. SUMMARY AND CONCLUSION 

In this work we have calculated, the neutrino mean free path with NFL corrections in 
quark matter for both degenerate and nondegenerate neutrinos. Contribution to the total 
mean free path arises from absorption and scattering processes. We find, for absorption 
processes, corrections to the MFP is O(tt^) for degenerate neutrinos and for non-degenerate 
case it is 0{al). This is due to the fact, that for non-degenerate neutrinos, quarks and 
electrons are coUinear in momenta on their Fermi surfaces, for which phase space for direct 
URCA processes vanishes. To get nonvanishing phase space which requires that the Fermi 
momenta of quarks and electrons should satisfy triangular relation, one must take into 
account the interaction between the quarks. 

Armed with the results of the previous sections, we now estimate the numerical values 
of the neutrino mean free paths in quark matter both for degenerate and nondegenerate 
neutrinos. The mean free paths are listed in Table-|T] in unit of meter where the neutrino 
energy {E,y) is set to be equal to the temperature (T) and rrig = 10 Mevj4|. Following ref.jst 
we take the quark chemical potential fiq ~ 500 MeV corresponding to densities pb ~ 6po, po 
is the nuclear matter saturation density. The electron chemical potential is determined by 
using the charge neutrality and beta equilibrium conditions which yields /ie = H Mev. The 
other parameter is used here are same as jsj. 

From Table-Hl we see that the anomalous logarithmic terms dominate the mean free path 
and reduces its value considerably. This reduced mean free path results in higher rate of 
scattering leading to higher rate of cooling of the compact stars for nondegenerate neutrinos. 
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We find the relatively short neutrino mean free path for degenerate neutrinos. These suggest 
that conditions are favorable for building up of the neutrino number in the core and the 
neutrino distribution function may approach close to that required for the enhancement of 
local thermodynamic equilibrium. 
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Recently, there has been a substantial effort to study the properties of cold and warm 
quark matter. Such studies are important to understand the properties of the astrophysical 
' compact objects like neutron stars and pulsars. There has been lot of experimental efforts 
^ \ like Einstein laboratory, ROSAT, CHANDRA and XMM, where various measurements are 
performed to understand the properties of neutron stars l|-|3| . 

There is a possibility that at the core of neutron stars the density may go upto 5 ~ 6 times 
the normal nuclear matter density where the matter is not expected to be in the hadronic 
phase 4, |5[. In fact, under such a scenario one expects that it would be more appropriate 
to describe the core of such dense stars as degenerate quark matter which is our main 
interest in the present work. 



It is known that the newly born neutron stars coo, 
antineutrinos within few minutes involving the direct 



Ui V J 



via the emissions of neutrinos and 
9| or the modified URCA processes 



[sl, [lo| . The direct URCA reactions can proceed in neutron rich matter if the ratio of the 
proton number density to the total baryon number density exceeds a critical value which 



2 



follows from the energy and momentum conservation properties In modified URCA 

process another nucleon catalyzes the reaction to occur under situations where the direct 
URCA reaction is forbidden. 

For quark matter, the dominant contribution to the emission of neutrinos is given by 
the quark analogus of /3 decay and the electron capture 9|]. These reactions are named 
as 'quark direct URCA' processes which have been studied in detail by Iwamoto 0, 
Our main focus here is to calculate the neutrino MFP, as mentioned earlier, for cold and 
warm QCD matter. Previously MFP of quark matter was derived in Ref.j?] where the 
calculations were restricted to the leading order by assuming free Fermi gas interactions. 
Here we plan to go beyond leading order to include 'plasma' or 'quasi particle' effects arising 
out of the interacting ground state. This, as we show, gives rise to logarithmic corrections 
;o the MFP in contrast to what one expects from the usual Fermi liquid theory (FLT) 
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131 ]. Similar non- Fermi liquid behavior of various thermodynamical quantities have 



recently drawn significant attention ^, s], [ol 1^ 16|. 



For example, in 



0, Q|; the authors have computed the leading contribution to the in- 
teraction part of specific heat and entropy when temperature is much smaller than the 
chemical potential of quark matter. It has been shown, in case of specific heat, for ideal 



gas the leading order term goes as (9(/i^T), while the correction term involve (9(T'^). In 
case of interacting CclSG, clS shown in ref . [lj| , the correction to the leading order contribution 
involve nonanalytic terms, known as non-Fermi liquid (NFL) corrections, which have been 



fi 



discussed extensively in j^, |5|, |9|, 



14-[l6|. Ref.[l5| examined NFL effects in the normal phase 
of high density QCD matter both using the Dyson-Schwinger equation and the renormaliza- 
tion group theory. Like specific heat, the magnetic susceptibility shows similar non-Fermi 
liquid behavior has been shown in 16j . Being motivated by these series of works we under- 
take the present investigation to estimate and see the consequences of such effects in the 
case of neutrino MFP. 



In FLT, quarks are treated as quas ipartic 



tional of the distribution function 



12 



13 



es and their energy (E) is regarded ClS Sb func- 



17H20|. FLT is restricted to the low-lying excita- 



tions near the Fermi surface, where the lifetime of quasiparticles are long enough. Therefore, 
it is an important tool to study the properties of nuclear (or quark) matter. In Ref. 2l|, it 
has been argued, that the exchange of dynamically screened transverse gluons introduces 
infrared divergences in the quark self-energies that lead to the breakdown of the Fermi liquid 



3 



description of cold and dense QCD in perturbation theory. A detailed study of non-Fermi 
liquid aspects of the normal state was presented in Ref.jsl. There the spectral density, 
dispersion relation and width of quasiparticles with momenta near the Fermi surface were 
derived at T = by implementing a renormalization group resummation of the leading 
logarithmic infrared divergences associated with the emission of soft dynamically screened 



transverse gluons 



We have already mentioned that such anomalous corrections are ultimately connected to 
the absence of magnetic screening of gluons via. Landau damping. One such calculations 
which we find to the most relevant for the present purpose was performed in j9|. It was 
shown that the emissivity receives logarithmic corrections is enhanced due to non-Fermi 
liquid effects. It might not be out of context here to recall two important works on neutrino 
mean free path in QED plasma. One is due to Tubbs and Schramm 2^ and the other is 



done by Lamb and Pethick |23|. In 22|, the resultant mean free path was calculated in the 
neutranized core and just outside the core. It is concluded in 23j, that neutrino degeneracy 
reduces the neutrino mean free path which suggest that neutrino may flow out of the core 
rather slowly. 

In our work we show that the neutrino mean free path receive logarithmic corrections 



where t 
the 



le dressed gluon propagator is used instead of bare propagator. Infact, we generalize 



23| results by incorporating the NFL corrections for the quark matter. The 



corrections to the MFP for degenerate and non-degenerate neutrinos, as we shall see, will 
involve different powers of 

In the interior of a neutron star, there are two distinct phenomenon for which the neutrino 
mean free path is calculated, one is absorption and the other one involves scattering of 



neutrinos \8 



a]. The corresponding mean free paths are denoted by l^ean ^mean- is to be 
noted that one also defines another MFP, known as transport mean free path which enters 
into the calculation of diffusion coefficient. The scattering MFP, on the other hand is related 
to the relaxation time that characterizes the rate of change of the neutrino distribution 
function |8|. One can combine /^"an "with to define total mean free path |24| . 

1 



Jtotal 
''mean 



1 1 

+ 



Jabs 
''mean 



Iscatt 
''mean 



(1) 



4 



II. MEAN FREE PATH 

In this section we calculate neutrino mean free paths in quark matter including NFL 
corrections both for degenerate and nondegenerate neutrinos. When the neutrino chemical 
potential (/ij/) is considered to be much larger than the temperature (T), the neutrinos 
become degenerate and for nondegenerate neutrinos fi^, <^ T. In our model the Lagrangian 
density is described by Ql 

Cw.ix) = ^l^{x)J^{x) + H.C. (2) 

where the weak coupling constant is G ~ 1.166 x 10"^^ in MeV units, and and are 
the lepton and hadron charged weak currents, respectively. The weak currents are 

l^{x) = 67^(1 -75)z/e + /i7^(l -75)1^^ + ...., (3) 
J^{x) = cos e,uY' {I -l5)d + sin e^uY - I5) s + (4) 



where 9^ is the Cabibbo angle (cos^ 9c — 0.948) 25|. 

The mean free path is determined by the quark neutrino interaction in dense quark matter 
via. weak processes. We consider the simplest (3 decay reactions; the absorption process 

d + z/g — 7- n + e~ (5) 

and other is its inverse relation 

M + — 7- (i + i/g (6) 

The neutrino mean free path is related to the total interaction rate due to neutrino 
emission averaged over the initial quark spins and summed over the final state phase space 



and spins. It is given by jsl 

1 g f d^pd 1 f d^pu 1 [ d^pe 1 



f d-'pd 1 r d'pu 1 [ d'% 1 .0 x4.4/p ,p_p_p\ 

ItLniEu.T) 2eJ (2vr)3 2Ej (2vr)3 2Ej (27r)3 2i5;e^^ ^ ^ ' " 

- n{pu)][l - n{pe)] - n{pu)n{pe)[l - n{pd)]}, (7) 



where, g is the spin and color degeneracy, which in the present case is considered to be 6. 
Here, E, p and Hp are the energy, momentum and distribution function for the corresponding 



5 



particle. |Mp is the squared invariant amplitude averaged over initial d quark spin and 
summed over final spins of u quark and electron as given by 



UG^cosH,{Pd-P,){P^-P, 



(8) 



Here, we work with the two flavor system as the interaction involving strange quark is 
Cabibbo suppressed [4, |9|. 



A. Degenerate Neutrinos 

We now consider the case of degenerate neutrinos i.e. when //^ ^ T or in other words we 
consider trapped neutrino matter. So in this case both the direct (Eqj5]) and inverse (Eqj6]) 
processes can occur and both the terms in Eq.([7]) under curly brackets are retained. 
Consequently, the (3 equilibrium condition becomes fid + fJ'u = /^w + /^e- Neglecting the 
quark-quark interactions and by using Eqs.([7]) and ([H]), for the mean free path one obtains 

= -^G^COS^ 6'c / d^Pd I d^Pu I d^Pe{l - cos 9du) (l - cos 9ue) 

X 6\Pd + P. - Pu - Pe)[l + e-^^''-^~^^^]n{pd)[l - n{p^)][l - n{p,)] (9) 

In the square bracket, the second term e'^^^"'^""^ is due to the inverse process (Eqj6]). Since 
the masses of u, d quark and electron are very small, one can neglect the mass effect on the 
mean free path. To carry out the momentum integration we define p = \pd + Pu\ = \Pu + Pe\ 



as a variable. Following the procedure, described by Iwamoto 



one has 



and 



pdp , . 

SYnBdudddv = — rjT — (10) 
Pf{d)pf{u) 

'l-cos9du){l-cos9ue) ^ . (\ (11) 

^Pi{d)pf[v)pf[u)pj[e) 



d^Pd = 2ir sin 9dud9duPldpd 

Pf{d) dpd 
= 2'K^--pdp-—dEd 
Pf{u) dEd 

= 2.^-I^W^du (12) 
Pf[v) duj 

d'p^ = 2nPl^^^dE,^du: (13) 
p doj 



6 



where we denote the single particle energy E^^u) as u. For the free case dp/du is the 
inverse quark velocity. It is well known that this slope of the dispersion relation changes in 
matter due to scattering from the Fermi surface and excitation of the Dirac vacuum. The 
modified dispersion relation can be obtained by computing the on-shell one-loop self-energy. 
For quasiparticles with momenta close to the Fermi momentum, the one-loop self-energy is 



dominated by the soft gluon exchanges 



26| . The quasiparticle energy cu satisfies the relation 



Ep{uj) + ReS(ci;,p(ci;)) 



(14) 



where we have approximated only to the real part of sel 



■ energy, since the imaginary part 



of S turns out to 
be obtained in 



De negligible compared to its real part 



27 



28 



Il6| . The detailed analysis can 



In the relativistic case, when the Fermi velocity f/ is close to the velocity of light c, 
exchange of magnetic gluons become important. In the non-relativistic case it is suppressed 
by a factor {v /cY with respect to exchange of electric gauge bosons and is usually neglected. 
The magnetic interaction, as stated before, is screened only dynamically and the problem 



can be written as 



14- 



21 



remains for the static gluons [26 
tion has to be performed 



26| 



29| . Therefore to obtain a finite result, a suitable resumma- 



. The analytical expressions for one-loop quark self-energy 



127r2 



iu - /i) In ( — ^ ) + i 

\u — u/ 



\UJ — jj\ 



(15) 

w — /i/ rZTT 

It exhibits a logarithmic singularity close to the Fermi surface i.e. when u ^ fi. Thus the 
on g ra nged character of the magnetic interactions spoils the normal Fermi-liquid behavior 



12 



13 



17 



20l |. The breakdown of the Fermi liquid picture is associated with the vanishing 
of the discontinuity of the distribution function at the Fermi surface js], [21|. This non- 
perturbative nature of the self energy gives rise to the non- Fermi liquid behavior. Here, m^) 
is a cut-off factor and should be an order of the Debye mass. Differentiating Eq. ffT^ with 
respect to p, we obtain at leading order in — is 



dp{(jj) 
dui 



1+^inr^^ 



37r 



p{u) 



(16) 



7 



where is the strong couphng constant, Cp = {N^ — l)/(2A'^c) and N^, is the color factor. 
Using Eqs. fll6p . (jH]) and flimi3p . the neutrino mean free path can be determined for two 
conditions. For \pf{u) —pf{e)\ > \pf{d) — pf{iy)\ 

1 



abs,D 
mean 



(17) 



To derive above Eq.f lT7|) . we use the resuh of integral js], Q 

roc rOG roo 

/ dEd / / dE,[l + e-^(^''-^'^V(P<i)[l - - niPe)]6{Ed + E, - E^ - E,) 

Jo Jo Jo 



Similarly, for \pf{d) — pf{iy)\ > \pf{u) — pf{e)\, the corresponding expression for mean 
free path can be obtained by replacing ^ f^d and /ig O //^ in Eq. lITTI) . Since quark and 
electron are assumed to be massless, the chemical equilibrium condition gives Pf{u)+Pf{e) = 
Pf{d) +pf{h'), which we use to derive Eq.l fTTl) . 

The other major contribution to the mean free path arises from quark-neutrino scattering. 
The neutrino scattering process from degenerate quarks are given by 



(19) 



for each quark component of flavor i{= u or d). The scattering mean free path of the 
neutrinos in degenerate case can be calculated similarly as evaluated by Lamb and Pethick 



m 



23| for electron- neutrino scattering. Assuming rriqjpf^ <^ 1 and including the non- Fermi 



liquid correction through phase space, the mean free path is given by 



1 



scatt,D 



I 



16 



ng^ao 



(E,-/i,)2 + 7r2T2 





[l + ^^"^ In 






I Sn 


K T )\ 



A(x, 



(20) 



and is the vector and axial vector coupling constant 
. In Eq.f l20p if we drop the color factor and the second 



Here, ruq- is the quark mass. C 
was given in TABLE-II of Ref. 
square bracketed term, we obtain results reported in 23| for dense and cold QED plasma 
with rriq replaced by rrie. In Eq. fl20|) the constants do = 4G^mg,/7r 22| and Uq- is the number 
density of quark, is given by 

d'^p 1 



n„ 



(27r)3e^(-^9«"^9.) 



(21) 



where 2 is the quark spin degeneracy factor. Exphcit form of A{xi) can be written as [8|, |23| 



4 Min(/it.,/igJ 

3 f^qi 



{C'+Ci){2 + -xi+2Cv^CA.x 



(22) 



and Xi = ^iu/Hq^ if ^-u < /^q, and Xi = figj^i^ if fiu > fJ'q, 



B. Nondegenerate Neutrinos 

We also derive mean free path for nondegenerate neutrinos i. e. when ^„ <^ T. For nonde- 
generate neutrinos the inverse process (jH]) is dropped. Hence we neglect the second term in 
the curly braces of Eq.(I7]). In this case, only those fermions whose momenta lie close to their 
respective Fermi surfaces can take part in a reaction. It is to be mentioned here, if quarks 



are treated as free, as discussed in 18 



33 



34j . the matrix element vanishes, since u, d quarks 
and electrons are collinear in momenta. The inclusion of strong interactions between quarks 
relaxes these kinematic restrictions resulting in a nonvanishing squared matrix amplitude. 
Since the neutrinos are produced thermally, we neglect the neutrino momentum in energy- 
momentum conservation relation js]. This is not the case for degenerate neutrinos where 
Py^T and therefore such approximation is not valid there. By doing angular average over 
;he direction of the outgoing neutrino, from Eq.(|8]) the squared matrix element is given by 
9| 



|M|2 = COS^ e,p){yd ■ Pu){Vu ■ Pe) 

= 6AG^ cosH,p]E,fi,^^ (23) 



where V = {l,Vf) is the four velocity. To calculate Eq. fl23|) we have used the chemical 
equilibrium condition fid = fJ'u + f^e and also the relations derived from Fermi liquid theory 



are given byjo] 



VF = l- — — ; 5n = fi (24) 

ZTT TT 



Putting |Mp in Eq.Q we have. 



1 ^^^""ic'cos' 9 J d'pd I d'p J d'p^ 



iabs,ND A-ttG 
•'mean 

5\Pd + P. - Pu - Pe)n{pd)[l - n{p.,)][l - n{p,)] (25) 
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Neglecting the neutrino momentum in the neutrino momentum conserving 5 function, the 



integrals can be decoupled into two parts. Following the procedure of Iwamoto 
angular integral is given by 



A 



the 



(26) 



and the other part 



B 



>dpd 



■dEd 



2 dpu j„ 



pldE, 



6{Ed + E,-Eu- Ee)n{pd)[l - n{pu)][l - n{p,)] 



(27) 



Changing the variables to Xd = {Ed — Hd)f^, Xu = —{E^ — jJiu)P and Xe = —{E^ — /ie)/3 
and denoting the single particle energy E^i^d) as u we have from Eg. ( 1271) 

dpd{uj) dpu{uj) 2 2 2; 



B 



duj du 



-PdPuPe^i^d + Xu + Xe + [3 E^)n{x d)n{-Xu)n{-Xe) (28) 



As the contribution dominates near the Fermi surfaces, extension of lower limit is a reason- 



able approximation 32 



Using Eq. f|T6|l and perform the integration of Eq. fl28l) following the procedure defined in 



16 



32 



35| . we have 



1 + ^lnf^) 
Sn \ T J 



2(1 + e-/5^' 

Using Eqs. (!25l) . (!26|) and (!29|l . the mean free path at leading order in T//i is given by 
1 



(29) 



: — G COS 9c l^d fi^u Aie Q _^ e-^^-) I- 



(30) 



TT- (l + e-z^-^j L 37r V T 7. 

The first term is known from js] and the additional terms are higher order corrections to 
the previous results derived in the present work. 

For the scattering of nondegenerate neutrinos in quark matter, the expression of mean free 



path was given by Iwamoto [8 



We incorporate the anomalous effect which enters through 



phase space modification giving rise to 



iscatt,ND on ^ \ m I \ ii 

Imean ■^^ ^'"'qi^ ^ H'l 

Here, we have assumed mg-fpf- ^ 1 and the constants ctq and number density rig. defined 
earlier. 



20 



(E,\^(E, 



CpOs 

37r 



In 



T ) 



(31) 
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III. RESULTS AND CONCLUSION 



Armed with the results of the previous sections, we now estimate the numerical values of 
the neutrino mean free paths. Here, E,y is set to be equal to 3T and ruq = 10 Mevjs], [2^. 
For the quark chemical potential, following ref. jsl, we take /ig ~ 500 MeV corresponding to 
densities pb ~ Qpo, where po is the nuclear matter saturation density. The electron chemical 
potential is determined by using the charge neutrality and beta equilibrium conditions which 
yields pe = 11 Mev. The other parameters used are same as j^. 



25 
20 



B 15- 
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with NFL correction 
without NFL correction 
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FIG. 1: Neutrino mean free path in quark matter for degenerate neutrinos. 



From fig ([T]), we find for degenerate neutrinos the anomalous logarithmic terms reduces 
the value of the mean free path appreciably both in the low and high temperature regime. Fig 
([2]) shows that for nondegenerate neutrinos NFL corrections is quite large at low temperature 
while at higher temperature it tends to merge. It is interesting to see from these two plots, 
that NFL corrections to the MFP in degenerate neutrinos is less than that of nondegenerate 
neutrinos. This reduced mean free path is expected to influence the cooling of the compact 
stars. It is also to be noted that all the results presented above are restricted to the leading 
log approximation which can be extended to take the next to leading order corrections into 



account 



36[. We plan to undertake such investigations in a future publication [37|. 
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FIG. 2: Neutrino mean free path in quark matter for nondegenerate neutrinos. 
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